On Hamiltonians Generating Optimal-Speed Evolutions 
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We present a simple derivation of the formula for the Hamiltonian operator(s) that achieve the 
fastest possible unitary evolution between given initial and final states. We discuss how this formula 
is modified in pseudo-Hermitian quantum mechanics and provide an explicit expression for the most 
general optimal-speed quasi-Hermitian Hamiltonian. Our approach allows for an explicit description 
of the metric- (inner product-) dependence of the lower bound on the travel time and the universality 
(metric-independence) of the upper bound on the speed of unitary evolutions. 
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In quantum mechanics the travel time for unitary time- 
evolutions between an initial and a final states A/ and A_f 
has a minimum that is proportional to the distance be- 
tween A/ and A_f in the state space JL] • For a system with 
a fixed energy scale, this implies that the speed of unitary 
evolutions has an upper bound. The problem of deter- 
mining a Hamiltonian operator that achieves the high- 
est evolution speed has been considered in [1, S 13, S ■ 
The purpose of the present article is to offer a very simple 
derivation of the formula for a time- independent optimal- 
speed Hamiltonian that can be directly generalized to the 
seudo-Hermitian representation of quantum mechanics 
In particular, we give, for the first time, the ex- 
plicit form of the most general time-independent quasi- 
Hermitian 8] optimal-speed Hamiltonian that evolves A/ 
into \p. 

In the standard formulation of quantum mechanics, the 
(pure) states of a physical system are identified with the 
rays in a complex separable Hilbert space Ti. Each ray 
can be determined in terms of an associated state vector 
tp GTi— {0} according to A^ — {ct/'|c G C}. It is usually 
convenient to use the one-to-one correspondence between 
the states A^ and the projection operators A^ := fli\L\ 
to describe the properties of the space of all states, i.e., 
the projective Hilbert space V{Ti.). For a Hilbert space 
Ti of dimension N < oo, 7^(7i) is the complex projective 
space CP^~^ that plays a central role in the description 
of geometric phases @ . 

It is an easy exercise to show that A^ satisfies 



A'=Av,=A]^, tr(A^) = l, 



(1) 



where "tr" denotes the trace. Recall that for a linear 
operator L acting in H, tr(L) :— J2n=i{^n\L ^n) , where 
{^n} is an arbitrary orthonormal basis of ?i. 20]. 

In view of ([1]), A^ is an element of the space *82(7i) 
of all linear operators L -.Ti ^ Ti fulfilling tr(L'f L) < oo. 
We can use "tr" to define the following inner product 
on this space {L\J) :— tr(L^^J). This makes 052 (7i) into 
a separable Hilbert space known as the Hilbert-Schmidt 
class [l3|- Because the state space V{T-L) is a subset of 
5B2(H), we can use the inner product (•!•) to define a 
notion of distance (metric) on P{'H). We define the line 



element ds on V{T-C) according to 



ds^ 



{dA,f,\dA^) 



{ib\ib){d^\d^p) - md^p)\' 



(2) 



where we have used A^ 



.V,:=™and(IT]), [m. ForiV< 
oo where ip can be represented by a complex column vec- 
tor 3* with components 3i,32: ■ • • ,3^, © takes the form 

^■5^ = Yl!l,b=l9ab'dladll where Qab' ■= \l\~'^{W^ab - 

■^*„}ih )- This is precisely the Fubini-Study metric tensor 
[12 |. For N = 2, endowing V{TC) with this metric yields 
a round two-dimensional sphere of unit diameter. 

Now, suppose that we wish to use an arbitrary Her- 
mitian (self-adjoint) Hamiltonian operator H : H ^ H 
to evolve an initial state X^^ to a final state A^^ . We 
can view the evolving state X-^m as a point moving on 
V{'H). According to ([2]), the instantaneous speed of the 
evolution is given by 



ds y/(i;it)\'4,{t)){^itMit)~ 
dt ~ MtWit)) 



liHtMim 



AE, 



i'it) 



where 



AE. 



^{t) 



'm)\H^^{t)) \m)\H^{t)w 



(^(i)|V.(i)) 



m)mtw 



(3) 



is the uncertainty in energy and we have employed the 
Schrodinger equation, H%l}{t) = itnp{t). We can integrate 
ds/ dt to obtain the length of the curve traced by X.^(t) in 
V{Ti) as a function of the travel time r, [i|: 



1 



^E.^{t) dt. 



(4) 



Because AE^(t) > for all t e [0, r), s is a monotonically 
increasing function of r. This makes r a monotonically 
increasing function of s. Therefore, the shortest travel 
time is achieved for the paths of the shortest length, i.e., 
the geodesies on V{Ti), [J]. 

Note that the geodesic distance is uniquely determined 
by the initial and final states and is insensitive to the 
choice of the Hamiltonian one uses to evolve the initial 



state along such a geodesic. The travel time depends on 
the Hamiltonian through the energy uncertainty AE^uy 
In particular, if one can make the latter arbitrarily large, 
the travel time can be made arbitrarily small. In typical 
situations, however, AE^u\ has a constant upper bound. 
For example, consider the case that the Hilbert space is 
finite-dimensional (N < oo) and the energy eigenvalues 
En are bounded functions of time; there is some £ G R+ 
such |£'„(i)| < £ for all n and t. Then, we can easily 
show that AiJ^(() < £; the travel speed is bounded by 
£/h; and the travel time has hs/£ as a lower bound. Here 
s is to be identified with the geodesic distance between 
the initial and final states. 

The above argument is valid, if one docs not have ad- 
ditional restrictions on the choice of the Hamiltonian. In 
practice, one may have to impose constraints that would 
make it impossible to evolve the initial state along the 
shortest geodesic connecting it to the final state. In this 
case one can formulate the problem as a constrained vari- 
ational problem [5|] . In the remainder of this article we 
consider constant unconstrained Hamiltonians where the 
minimum travel time depends, besides the geodesic dis- 
tance between the initial and final states, on a single real 
parameter specifying the energy scale of the system. 

Let H he a, time-independent Hamiltonian operator. 
Then the time-evolution operator e~'*^/^ commutes with 
H and H^, and A£'^(t) does not depend on t. In this case, 
(|4]) implies r — hs/AE^, and the speed of the evolution 
is given by AE^/h. Therefore, to achieve the highest 
speed we need to choose the Hamiltonian so that AE^/h 
is maximized. This shows that the travel time is bounded 
by the ratio of the minimum of s, i.e., the geodesic dis- 
tance between A^^ and A^^, to the maximum of speed 
AE^/h. 

Because we require the evolving state A^(4) to trace 
a geodesic in V(H) that connects A^^ and A^^^, it lies 
entirely in the projective Hilbert space V{'H') where H' 
is the subspace of H spanned by ipj and ipp ■ This is in 
fact a characteristic property of the Fubini-Study metric 
[l|, [l^l ■ It shows that we can restrict our attention to the 
case that TC is two-dimensional; N = 2, [J]. Furthermore, 
without loss of generality, we can suppose that tT{H) — 0. 
This implies that the eigenvalues of H have opposite sign, 
E2 — —El =: E. Let {V'ij''/'2} be an orthonormal basis 
consisting of the eigenvectors of H, Hipn — Enipn- We 
expand "0(0) = ipi in this basis to find 



i^I = CiV'l + C2V'2 



ci,c2 e 



and use the time-independence of AE^ to compute it at 
i = 0. In view of ^ and (O, this yields 



AE, 



"0 



\CI?-\C2V 
|C1P + |C2|2 



< E. 



Therefore, the travel time r satisfies 

hs 



T>T„ 



E' 



(6) 



(7) 



where s is the geodesic distance between A^j^ and A^j^ 
in V{Tl). dl]) identifies Tmin with a lower bound on the 
travel time. Next, we construct a Hamiltonian H^, with 
eigenvalues ±E for which t — Tmin. This shows that 
indeed Tmin is the minimum travel time. 

Because s is completely determined by A.^^ and A^^, 
the condition r = Tmin is fulfilled if and only if AE^ = E. 
In light of ^ this is equivalent to |ci| — |c2|. If we 
expand ipp in the basis {f/'ii ^'2} to find 



tpF = dii}i+ (12^1^2, di,d2eC, 



(8) 



and compute AE^p at t — t, we obtain Q with (ci,C2) 
replaced with {di,d2). As a result, in order to maintain 
AE^p = E, we must have |di| = \d2\- 

Next we express |ci| — \c2\ and |di| = \d2\ in the form 
C2 = e'"jci and ^2 = e'^rdi, for some 0^,0^ G M, re- 
spectively. Substituting these in (O and ([5]), we find 

V>i+e*"^V2 = crV/, Vi+e*"^V'2 = drVf- (9) 

We can solve these equations for ipi and 7/12 in terms of 
ipi and ipF, and use the spectral resolution of iJ^, i.e.. 



iJ. =i?(-|V'l)(^l| + |^2)(V'2|), 



(10) 



to compute H^,. This calculation is more conveniently 
performed in terms of 



■d := a, —a. 



tpj 



tpi 



ijjp ■■= ^ , . (11) 



V2ci' ' V2di 

The result is [I, i] 

iE{\i^F){M-\^l){^F\) 



H, 



4sin(f) 
zgcot(f) ni;F){^i\ \i>i){M 

4 V {i'Ai^p) (V'fIV'/) 



(12) 



(13) 



where we have used the fact that 0/ and ipp £^re unit vec- 
tors. Moreover, m implies cos2(i?/2) - (^Jfgt^'f^,) , 



which as explained in [l| identifies d with 2s, [2l| 

Equation p^ can be easily modified to give the ex- 
pression for the optimal-speed Hamiltonians in pseudo- 
Hermitian quantum mechanics. One merely needs to 
make the following substitution in the above analysis 



(5) IV-n) -^ IV^n^, (V'nl ^ ^ "^n | := (V""!??-!-, 



'i)+i 



(14) 



where 77+ : 7i — > 7i is the metric operator that defines 
the inner product of the physical Hilbert space Tiphysj i-e., 
(•, •)^_|_ :— (•|77+-) —^ 'l'^; a-nd s,,^ is the distance defined 
by the natural metric on the projective Hilbert space 
■p(7ipiiys)- We can obtain the line element associated 
with this metric by making the substitutions (fT4|l on the 
right-hand side of ([2]). This gives [ll| 



ds 



v+ 



^'ijj\'ijj>-^d'ip\dip> 1 ^VMV'^ 



(15) 



Similarly we find the following expressions for mini- this relation in ((T5| and using ([M]) and (fT9l) we find [Tl| 



mum travel time Tj^^jn and the optimal-speed ry+-pseudo- 



Hermitian \i£ 

±E). 



Hamiltonian H^'^ (with eigenvalues 



^-7+) 



hs 



v+ 



E 



m" 






where 



cos^(s,,^) 



^V/IV'-F^ 



^'i(ji\ipiy^'il)F\^F'^' 



(16) 
,(17) 



(18) 



According to ([7]) and (fT6)) . if wc choose 77+ such that the 
geodesic distance s,,_^ between A^^ and A^^^ in 7^(7iphys) 
is smaller than the geodesic distance s between A^^ and 

X^p in V{Ti), then t^^ < Tmin. This is the essence of 
the main result of [Q- Indeed, as we show below, it 
is possible to choose 77+ so that regardless of the choice 
of A^j and A^^ their distance in 7^(7iphys) becomes arbi- 
trarily small. But this does not seem to have any physical 
implications, for such an evolution amounts to evolving 
a state to an arbitrarily close state in an arbitrarily short 
time. The physical quantity of practical significance, par- 
ticularly in areas such as quantum computation, is the 
speed of the evolution, namely E/h, which is a univer- 
sal quantity independent of the choice of 77+. Therefore, 
the minimum travel time between states of a given dis- 
tance is independent of 77+, [11]. A physical process that 
involves evolving A^,^ into A^,^ in 7^(7iphys) using an ?/+- 
pseudo-Hermitian Hamiltonian iJ : 7i — > 7i in time r 
may be described equally well by evolving A 1/2 , into 



A 1 



±/2 , 



in V(Ti.) using the equivalent Hcrmitian Hamil- 



1/2 



tonian h := ri^ H r]^ in the same time r. As shown 
in [11|, the length of the curve corresponding to these 
evolutions in the respective projective Hilbert spaces are 
identical. Therefore, they will have the same speed. 

Next, we wish to show how by choosing the metric 
operator we may adjust the value of s,,_|_ and consequently 



.(^7+) 



Again, without loss of generality we confine our 



attention to the case N 
standard basis of C^ 



2. Let {61,62} denote the 



I.e., 61 := 



62 



Then 



we can represent any metric operator t^^ in {61, 62} by a 
positive-definite matrix of the form 



V+ 



a b* 
b c 



(19) 



where a, c G K and fa e C. Because 7y+ is a positive- 
definite matrix, 

a+c = tr{rj+) > 0, D := ac-\bf = det(^) > 0. (20) 

For states A^ , differing from Aej , we can use a representa- 
tive state vector of the form ip :—( }^ . ] . Substituting 



ds 



Didx"^ +dy'^ 



')+ 



[a + 2{bix + 622/) + c(.t2 + 2/2)]- 



(21) 



where 61 and 62 are respectively the real and imaginary 
parts of fa, i.e., fa =: 61 -I- ib2. 

Next, we introduce the angular coordinates {(p, 9) that 
are related to {x, y) according to x = tan(|) cos((p + (3), 
y = tan(|)sin((^ -I- j3), where (3 := tan~^(62/fei)- These 
coordinates also allow for treating the state Aga- To see 
this, first observe that for {x, y) € M^, we have ip G [0, 27r) 
and 9 £ [0,7r). The state X^^ corresponds to the point 
at infinity in the x-y plane which we can identify with 
9 = TT. In terms of {ip, 9), (|2ip reads 



ds 



fci {d9^ 



sm 



'V 



v+ 



[1 + fc2 cos 9 + k^ cos ip sin 9] ^ 



(22) 



D 



tr(j)+)^ ' ''^2 



where we have introduced fci 

^-p^, and k-i := — L-L. Note that because of (1^01) we have 

a+c ' ^ a+c ' r 

fci > 0, -1 < /c2 < 1 and < k^ < 1. 

If we set 7]+ = I we recover the standard Euclidean 
inner product on the Hilbert space. In this case 3^+ = s, 
a == c = 1, fa = 0, fci = 1/4, k2 = k-i = 0, and ^ 
becomes ds^ — j{d9'^ + sir? 9 dp^) . This is just the stan- 
dard metric for a round sphere of unit diameter. In |ll| 
we show that 7'(7iphys) is related to ViTi) by an isome- 
try. Therefore, (|22p also describes a round sphere of unit 
diameter, and [p, 9) are the usual spherical coordinates. 

We can use ([T7| and p^ . to obtain the explicit form 
of the optimal-speed 77+-pseudo-Hermitian Hamiltonians 
for given initial A^j^ and final A^^ states. We can always 
perform an invertible linear (basis) transformation in Ti 



so that A 



Tpl 



A. 



Let X^p be an arbitrary final state 



(that is different from Ae^ ) . Then we can take 7/;/ 



and tpi 



for some C G C In view of P^ and ([T 



we have <^i\il)i)~= a, ~i'i(jp\'iljpy= — ^^i-, -^ipjlt/jpy— 
aC -I- fa* =: 1^. Inserting these in (|18p and using P^ . we 
find coss^+ = \^\/\/dTW and 



Jv+) 



-^ cos-i (l^l/v/^T^ 



(23) 



Similarly we employ (flTl) and p9|) . to obtain the ma- 
trix representation of HiP^ in the basis {61, 62}: 



i/^('?+) = 



iEe- 



4a^/D 



-afa* 

„2 



ab* 



,(24) 



where u; :— arg(^), i.e., 6^^ — £,/\£,\. It is interesting 
to see that the minimum travel time and optimal-speed 
Hamiltonians are respectively determined by the modu- 
lus and the phase of ^. Note also that the right-hand side 
of ([M)) does not have a unique limit as ^ — > 0. This is 
because ^ = corresponds to the case that A^^ and A^^, 



are antipodal points of PCHphys) that are connected via 
an infinity of geodesies with equal length. 

Setting a = c = Z) = 1 and b == in ^ and ^, we 
find the explicit form of the optimal-speed Hamiltonian 
and the minimum travel time in conventional quantum 

mechanics: H±= ^ i ^-iu, q 

^min - I COS-1 (iCI/yiT^) . (25) 

Note that in this case ^ = C ^-nd e*" = C/ICI- Com- 
paring ([25|) with ([23]), we see that by keeping a and b 
fixed, so that ^ is left unchanged, and decreasing the 
value of c we can make £) as small as we wish. This in 



turn reduces the value of r. 



(»)+ 



below that of r„ 



For 



example, we can set a = \ and b 
^ = C, and we find iJ*(''+) = 



4 




and 



J'7+) 



I cos-1 



ni^^^w 



For c < 1, this 



J'?+) 



yields t^^ < Tmin- As we explained above, this ob- 
servation does not seem to have any practical implica- 
tions, if we use Hi'^ to generate a unitary evolution, 
i.e., consider the dynamics taking place in ViHphys)- If 

we instead consider the dynamics defined by -ff*'^ in 
V{'H), then the travel time is still given by (|^5|) (which 
can be made smaller than T,„in) but the evolution is non- 
unitary. This is a manifestation of the nonexistence of 
an upper bound on evolution speed for non-unitary evo- 
lutions [l6|. A more interesting observation is that one 
can realize this fact using a quasi-Hermitian Hamilto- 
nian (T'T-symmetric Hamiltonians considered in 15] be- 
ing special cases) (171 : there is no upper bound on the 
speed of quasi-unitary evolutions, [22 • 

To offer a physical interpretation for this result we first 



recall that in quantum mechanics, a physical system is 
represented by a Hilbert space-Hamiltonian pair (7i, H). 
This representation is however not unique, for unitary- 
equivalent Hilbert space-Hamiltonian pairs describe the 
same system. If Ti and Tiphys are Hilbert spaces that 
have identical vector space structure but different inner 
products (say corresponding to the choices / and ?]+ for 
their metric operators respectively) , one can use a single 
Hamiltonian operator to represent two different quantum 



Tiri+) 



('l+)^ 



systems, e.g., {K,Hi^') and (7^phys,-ff* ) represent 
distinct physical systems with quasi-unitary and unitary 
dynamical evolutions, respectively. The above argument 
shows that while the evolution speed for the latter sys- 
tem is given by E /h, that of the former can be made 
arbitrarily large. To determine whether this observation 
can have practical applications requires a more detailed 
investigation of the role of quasi- and pseudo-Hermitian 
Hamiltonians in open quantum systems [l8|. 

In summary, we have offered a straightforward deriva- 
tion of an expression for the most general time- 
independent optimal-speed quasi-Hermitian (in particu- 
lar Hermitian) Hamiltonians and established by explicit 
calculation the metric-dependence of the minimum travel 
time and metric-independence of the maximum travel 
speed. Our analysis confirms the existence of infinitely 
fast quasi-unitary evolutions. These might find applica- 
tions in areas such as quantum computation and quan- 
tum control, [3|, [l9| . The derivation of an explicit expres- 
sion for the most general optimal-speed quasi-Hermitian 
Hamiltonian, that we have reported here, is a necessary 
step in this direction. 

Acknowledgments: This project was supported by the 
Scientific and Technological Research Council of Turkey 
(TUBITAK) as a part of the 2007 TUBITAK Science 
Award in Basic Sciences, and by the Turkish Academy 
of Sciences (TUBA). 



[1] J. Anandan and Y. Aharonov, Phys. Rev. Lett. 65, 1697 

(1990). [10] 

[2] G. N. Fleming, Nuovo Cimento A 16, 232 (1973); 

L. Vaidman, Amer. J. Phys. 60, 182 (1992). [11] 

[3] N. Margolus and L. B. Levitin, Physica D 120, 188 [12] 

(1998). 
[4] D. C. Brody, J. Phys. A 36, 5587 (2003). [13] 

[5] A. Carlini, A. Hosoya, T. Koike, and Y. Okudaira, Phys. 

Rev. Lett. 96, 060503 (2006); Phys. Rev. A 75, 042308 [14] 

(2007). 
[6] D. C. Brody and D. W. Hook, J. Phys. A 39, L167 (2006). [15] 

[7] A. Mostafazadeh and A. Batal, J. Phys. A 37, 11645 

(2004). [16] 

[8] F. G. Scholtz, H. B. Geyer, and F. J. W. Hahne, Ann. 

Phys. (NY) 213 74 (1992). [17] 

[9] D. N. Page, Phys. Rev. A 36, 3479 (1987); A. Bohm [18] 

and A. Mostafazadeh, J. Math. Phys. 35, 1463 (1994); 

A. Bohm, A. Mostafazadeh, H. Koizumi, Q. Niu, and [19] 

J. Zwanziger, The Geometric Phase in Quantum Sys- 



tem,s, Springer, New York, 2003. 

M. Reed and B. Simon, Functional Analysis, vol. \, Aca- 
demic Press, San Diego, 1980. 

A. Mostafazadeh, Phys. Rev. Lett. 99, 130502 (2007). 
T. Eguchi, P. B. Gilkey, and A. J. Hanson, Phys. Rep. 66, 
213 (1980). 

I. Bengtsson and K. Zyczkowski, Geometry of Quantum 
States, Cambridge University Press, Cambridge, 2006. 
A. Mostafazadeh, J. Math. Phys. 43, 205, 2814, and 3944 
(2002). 

C. M. Bender, D. C. Brody, H. F. Jones, and B. K. Meis- 
ter, Phys. Rev. Lett. 98, 040403 (2007). 
P. E. G. Assis and A. Fring, J. Phys . A 41, 244002 (2008); 
A. Mostafazadeh, arXiv:07 09l756] 
H. F. Jones, private communications. 
M. Jakob and S. Stenholm, Phys. Rev. A 70, 012104 
(2004). 

N. Khaneja, R. Brocket!, and S. J. Glaser, Phys. Rev. A 
63, 032308 (2001). 



[20] tr(L) is independent of the choice of {^,i}. 
[21] In view of this expression for cos^(i9/2), (|lip . and 

= 1, ci = V(V'/i^/>/2e'^ and 
where 7, 5 G R are arbitrary if 



{'IpFli'F} 



di = V(^i^lV'F>/2e'* 



(i/)i|i/'f) = 0, otherwise e'^*"^' = eV(^^|^j,)/|(^^[^^)|. 
[22] U : Ti. ^ TC is quasi-unitary, if it is ?;+-pseudo-unitary 
for some metric operator ri+, i.e., U~^ = ri+U^r]^^. See 
A. Mostafazadeh, J. Math. Phys., 45, 932 (2004). 



